%------------------Ormsby子波的时域表达式---------------%
function w = Ormsby(dt,N,fL1,fL2,fH1,fH2,fa)
    t=-N/2*dt:dt:N/2*dt;
    s1= (((pi*fH2.*sinc(pi*fH2*t)).^2 - (pi*fH1.*sinc(pi*fH1*t)).^2)/(pi*fH2-pi*fH1))...
            - (((pi*fL2.*sinc(pi*fL2*t)).^2 - (pi*fL1.*sinc(pi*fL1*t)).^2)/(pi*fL2-pi*fL1));
    s2=hilbert(s1);
    w=real(s2)*cos(fa)+imag(s2)*sin(fa);
end